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Summary! For general relativistic spacetimes filled with an irrotational perfect 
fluid a generalized form of Friedmann's equations governing the expansion factor of spa- 
tially averaged portions of inhomogeneous cosmologies is derived. The averaging prob- 
lem for scalar quantities is condensed into the problem of finding an 'effective equation 
of state' including kinematical as well as dynamical 'backreaction' terms that measure 
the departure from a standard FLRW cosmology. Applications of the averaged models 
are outlined including radiation-dominated and scalar field cosmologies (inflationary 
and dilaton/string cosmologies). In particular, the averaged equations show that the 
averaged scalar curvature must generically change in the course of structure formation, 
that an averaged inhomogeneous radiation cosmos does not follow the evolution of the 
standard homogeneous-isotropic model, and that an averaged inhomogeneous perfect 
fluid features kinematical 'backreaction' terms that, in some cases, act like a free scalar 
field source. The free scalar field (dilaton) itself, modelled by a 'stiff' fluid, is singled out 
as a special inhomogeneous case where the averaged equations assume a simple form. 



1. Introduction 



The present paper continues a line of research on average properties of inhomogeneous flu- 
ids in general relativity that is based on a simple and intuitive averaging procedure. The 
simplification is guided by the restriction to scalar dynamical variables and to standard 
volume integration, in which case averaging is straightforward (for a discussion of alterna- 
tive procedures see, e.g., Stoeger et al. 1999). Averaging is aimed at the construction of an 
effective dynamics of spatial portions of the Universe from which, in principle, observable 
average characteristics can be inferred like Hubble's constant, the effective 3-Ricci scalar 
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curvature and the mean density of a given spatial domain, which is bounded by the limits of 
observation. Naturally, this view entails a scale-dependent description of inhomogeneous 
cosmologies. In the case where the extension of the (simply-connected) spatial domain to 
the whole Universe is possible, such a description may allow to draw conclusions about 
global properties of the world models. 

Paper I (Buchert 2000) was concerned with 'dust cosmologies' restricting attention 
to the most popular inhomogeneous cosmologies. It is, however, desirable to extend the 
regime of application of an 'effective' (i.e. averaged) dynamics to a wider range of spatial 
and temporal scales than that covered by the matter model 'dust'. This is the motivation 
of the present work which presents the results for a large class of perfect fluid cosmologies. 

This class opens quite a piece of new terrain: it covers radiation-dominated cos- 
mologies, scalar field cosmologies including inhomogeneous dilaton/string cosmologies and 
inflationary cosmologies. It also extends the range of validity concerning averages of large- 
scale structure formation models for collisionless matter, in which case the presence of a 
pressure-force that counteracts gravity is implied by the development of multi-streaming 
within high-density regions; here, it provides a phenomenological extension by including 
physics on smaller spatial scales for the evolution of structure (see: Buchert & Dommguez 
1998, Buchert et al. 1999 in Newtonian cosmology; Maartens et al. 1999 in GR). 

This paper is organized as follows. Sect. 2 presents Einstein's equations for irrotational 
perfect fluids with the choice of foliation into flow-orthogonal hyper surfaces. Averaging 
the scalar parts of Einstein's equations is investigated in Sect. 3. The result is presented 
in a Theorem in Subsect. 3.2., which shows that the average expansion of inhomogeneous 
models is controlled by 'kinematical backreaction' due to shear and expansion fluctuations, 
and by 'dynamical backreaction' due to a non-vanishing pressure gradient in the hypersur- 
faces. It is manifest that the simple relation between averaged 3-Ricci scalar curvature and 
'kinematical backreaction', as found for the matter model 'dust' in Paper I (Buchert 2000), 
is supplemented by several additional effects: besides dynamical contributions to 'backre- 
action' the averaged energy- and momentum conservation laws do not yield conservation 
laws for the averaged fields. This leads to more drastic changes in the average flow com- 
pared with the standard model. Corollary 1 and Corollary 2 present compact formulations 
of the averaged equations for effective sources. The relations between additional sources in 
the generalized Friedmann equations are so formally reduced to the search for an effective 
equation of state. While Sect. 3 appears rather formal, especially because the presence 
of pressure involves an inhomogeneous lapse function and so impairs the simplicity of the 
equations, emphasis is focused on the application side in Sect. 4. There we discuss some 
relevant subcases that are members of the family of barotropic fluids: averaged 'dust' mod- 
els are recovered from the more general framework, averaged radiation-dominated models 
display deviations from a standard radiation cosmos, even if 'kinematical backreaction' is 
absent, and the application to inhomogeneous scalar field cosmologies is outlined. 
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2. Einstein's Equations for Perfect Fluids 



2.1. Choice of Foliation and Dynamical Variables 

We shall assume for the cosmic fluid that it is perfect and irrotational, so that we can 
introduce a foliation of spacetime into hypersurfaces orthogonal to the 4-velocity. It is 
not a problem to allow for a 'tilted' slicing in order to include, e.g., vorticity (see, e.g., 
MacCallum & Taub (1972), King & Ellis (1973), Hwang & Vishniac (1990), and Ellis et al. 
(1990)). For the applications we have in mind and also to keep the present investigation 
transparent, we shall evaluate everything for this class of fluids. 

For the purpose of averaging we shall consider a compact and simply-connected do- 
main contained within spatial hypersurfaces that are specified below. This domain will be 
followed along the flow lines of the fluid elements; thus we require that the total restmass 
of the fluid within the domain be conserved. 

Let us first consider the (conserved) restmass flux vector 1 

ikP := gu» ; = ; g > , (la) 

where g is the restmass density and the flow lines are integral curves of the 4-velocity -u M . 
Confining ourselves to irrotational fluids guarantees the existence of a scalar function S, 
such that 

-d^S 

«"==—. 0»> 

where the function h will be identified below. It normalizes the 4-gradient d^S so that 
u^u^ = -1, 



h = ^-d a Sd a S = u^d^S =:S>0 . (lc) 

The overdot stands for the material derivative operator along the flow lines of any tensor 
field T as defined covariantly by 

t:=u»T,» . (Id) 



We shall aim at a covariant description of the fluid flow with respect to the natural 
foliation of spacetime into hypersurfaces S = const, representing the 3-dimensional 'wave 
fronts' (for the covariant fluid approach compare Ellis & Bruni 1989, Bruni et al. 1990a,b, 
Dunsby et al. 1992). With our choice of the fluid's 4-velocity (lb) we have to assure that 
it remains time-like and, hence, the hypersurfaces S = const, space-like. For this to be 
true the 4-gradient of the scalar field has to be time-like, 

d a Sd a S = -h 2 < . (la) 



Greek indices run through 0...3, while latin indices run through 1...3; summation over 
repeated indices is understood. A semicolon will denote covariant derivative with respect 
to the 4-metric with signature (— , +, +, +); the units are such that c = 1. 
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(For h G IR this is always true.) As already noted, the definition (lb) implies that -u M is 
irrotational (2a) ; it also implies that the covariant spatial gradient of S in the hypersurfaces 
of constant S, denoted by S\\^, vanishes (2b), 

"V = h » ah fu[a;(3] = -h^hfi^d^S).^ = ; (2a) 

S\\p = h»d v S = d»S + u ll S = , (26) 

where := g^ v + u^u v is the projection tensor into the hypersurfaces S = const. 
orthogonal to the integral curves of the 4— velocity -u M , h^ v u u = 0. 

Hence, S(t) and h(X l ,t) play the role of 'phase' and 'amplitude' of the fluid's 'wave 
fronts', respectively. 

On these hypersurfaces we introduce the 3-metric gij (the first fundamental form) that is 
induced by the projection, as well as the extrinsic curvature tensor (the second fundamental 
form) : 

hij ■■= g tl uh^ i h v j = gij ; K tj := -u^ v h^ i h v J = -u i;j . (3a, 6) 

The final result will be covariant with respect to the given foliation, but we shall label 
the flow lines by introducing (intrinsic) Gaussian coordinates X 1 that appear in the line- 
element 

ds 2 = -N 2 dt 2 + g ij dX i dX j . (4a) 

Since by this choice of coordinates the velocities in 3-space vanish we are entitled to 
call X % Lagrangian coordinates. In the language of the ADM formalism, which is put 
into perspective in the Appendix, we have a vanishing shift vector in the hypersurfaces 
and the lapse function (together with the 3-metric) encodes the inhomogeneities. For 
scalar functions T = i/j the covariant derivative (Id) reduces to the total (or Lagrangian) 
derivative along the flow lines, 

d dr^ 1 

-^:=—d^ = u»d^ = -d t ^ , (46) 

where is the (inhomogeneous) lapse function. It is crucial to note that the latter operator 
corresponds to a total time derivative with respect to proper time r, which can be defined 
by r := / Ndt. 

For later discussions we may express the symmetric tensor K^, or the expansion ten- 
sor 0y := —K^, respectively, in terms of kinematical quantities and their scalar in- 
variants (Ehlers 1961). We decompose into its trace-free symmetric 'shear tensor' 
&ij '■= CTfiuh^^j , a^ v u u = 0, and its trace, the 'rate of expansion' 9 := u a . a . From the 
decomposition u^ u = a^ u + ^9h^ u — u^u^ we have: 

-K^ = <dij = a tj + ^9gij ; —K = 9 . (5a) 
The tensor has three principal scalar invariants; in what follows we shall use two of them: 

I := —K = 9 ; 211 := K 2 - K\K\ = -9 2 - 2a 2 , (56, c) 

3 
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where we have introduced the 'rate of shear' a by a 2 := \o % j(j J \. 



2.2. Basic Equations in 3 + 1 Form 

Einstein's equations for an irrotational perfect fluid with the energy-momentum tensor 

:= eupu v + phpy , (6) 

with energy density s and pressure p, may be cast into a set of 'constraint equations', the 
Hamiltonian and momentum constraints 2 , 

K + K 2 - K l 3 K\ = 16nGe , (7a) 

K i j\\i- K \i = ° ' W 
and 'evolution equations' for the the two fundamental forms: 

■j^-gij = -2 g ik K k j , (7c) 
-j^lCj = KK i j + n i j -A7rG8 i j (€-p)-(a i ] y + a i a j ) , {Id) 

where the acceleration is completely contained in the hypersurfaces of constant S and is 
defined as 

(H = h^a^ , a M := uV ;1/ = ii^ , a^u^ = , (8a) 

and 1Z := TZ l i , K := K % i denote the traces of the spatial Ricci tensor TZ l ■ and the extrinsic 
curvature tensor, respectively. Below, we shall only average the 4— divergence A of the 
acceleration field: 

A := a M ;At = a i ^ i + a l ai . (86) 

A nonvanishing acceleration is a consequence of the fact that the pressure term forces 
deviations from a geodesic flow. From T^ u v = we derive the energy and momentum 
conservation laws: 

u M 3% = e = -6(e + p) , (9a) 

a Ma T^ = & u a =a a = 1 —d ll ph» a = ^p, Q . (96) 

e + p e + p 1 



Hence 



1 . , 

ai = ■ — p\i , (9c) 

e + p 1 



2 As before, a double vertical slash abbreviates the covariant derivative with respect to 
the 3-metric g^; for scalars it reduces to the partial derivative with respect to Lagrangian 
coordinates denoted by a single vertical slash. 
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and 

From Eq. (la) we also have the continuity equation for the restmass density 

g + 9g = . (9e) 

According to (lc), S = -^d t S(t) = h, we can write Eq. (9d) completely in terms of the 
magnitude h and its spatial derivatives: 



A=[ ir)„r-h h " i+ v hh, ' = h {h ] ■ {10a) 

Two other derived formulas will be used in what follows. First, Raychaudhuri's equation, 
which follows by taking the trace of (7d) and inserting (7a): 

9 = -\e 2 -2a 2 -AnG{e + ?>p)+A , (106) 

and, second, an expression for the spatial Ricci curvature scalar in terms of the energy 
source terms, the restmass density, the magnitude h and its spatial derivatives: eliminating 
211 = §# 2 - 2a 2 from Eq. (10b) and, using Eq. (5c), from the Hamiltonian constraint Eq. 
(7a), we obtain with Eq. (9e): 

K =12. G ( £ - P )-4(I) +ft (I)'' . (10c) 



2.3. Thermodynamics of the Fluid 

First we note that the energy and restmass conservation laws Eqs. (9a,e) are equivalent 
according to the first law of thermodynamics, 

dp de ds , ^ 

= —— =■ — , 11a 

g e + p s 

upon dividing by dr. The latter equality defines the entropy density s that obeys the 
conservation law 

( S ^). M = ; s + 6s = . (lib) 

For closing the system of Einstein equations we need to identify a concrete matter model. 
Specific models are obtained by invoking a local 'equation of state' relating the pressure 
with the other dynamical variables. We shall discuss examples that are all members of the 
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class of 'barotropic fluids', i.e., p = a(e) is assumed to be locally given and, in particular, 
the function a is the same for each fluid element (at each trajectory). The special inho- 
mogeneous fluid cases discussed will all be contained in the simpler class a(e) = 7£ with 
7 = const., a 'dust' matter model (7 = 0), a 'radiation fluid' (7 = |), and a 'stiff fluid' 
corresponding to a free minimally coupled scalar field (7 = 1). 

We shall now identify the normalization amplitude h. Let us first derive h for a 
'barotropic fluid'. The momentum conservation law implies 



N h e+p 



(lib) 



Defining II := J with e = a~ 1 (p), w e may write Eq. (lib) as (h hi(^) + 1% = 0, 
which may be integrated to give 

h oc exp IT , (He) 



up to a time-dependent integration function. Hence, we may write 

dh dp 
h e + p 



(lid) 



In general, we identify the magnitude h with the 'injection energy per fluid element and 
unit restmass' (Israel 1976), 

h:= £ -±* , (He) 
Q 

which is related to the relativistic enthalpy r\ := by h = rj/m with m the unit restmass 
of a fluid element, and n the baryon density. Note that Eq. (lid) holds by defining h as 
in Eq. (lie) as a result of the conservation laws Eqs. (9a,e), since from Eq. (11a), 

de = hdg . (11/) 

For a barotropic fluid we can easily see that £ is a function of the restmass density only 
and, hence, h is a function of g. The evolution equation for h in this case (and in the 
simpler case p = 75) reads: 

h + ea'(e)h = h + 1 6h = $ ; (llg) 
h obeys a simple continuity equation in the case of a 'stiff' fluid with 7 = 1. 
The discussion of special cases will be resumed in Sect. 4. 
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3. The Averaged System 



3.1. The Averaging Procedure 

Spatially averaging equations for scalar fields is a covariant operation given a foliation of 
spacetime. Therefore, we shall in what follows only consider scalar functions T(X l , t). We 
shall define the averaging operation by the usual spatial volume average performed on an 
arbitrary compact support of the fluid V contained within the hypersurfaces S(t) = const.: 

(y) v := -L J^y Jd 3 X , J := ^det( 9ij ) . (12a) 

The volume of the region itself (set T = 1) is given by Vr>(t) := j v Jd 3 X. 

We also introduce a dimensionless scale factor via the volume (normalized by the volume 

of the initial domain Vr> ): 

<*(«) := (^) . (1») 

This means that we are only interested in the effective dynamics of the domain; ax> will 
be a functional of the domain's shape (dictated by the metric) and position. Since the 
domains follow the flow lines, the total restmass Mx> := f v gJd 3 X contained in a given 
domain is conserved. 

The following formulas are crucial for evaluating averages. Taking the trace of Eq. (7c), 
written in the form 

1 rl 

we obtain with \g lk -^9ki = (In J)' the identity 

j = 6J . (12c) 

The rate of change of the volume V(t) in the hypersurfaces S(t) = const, is evaluated by 
taking the partial time derivative of the volume and dividing by the volume. Since c\ and 
d 3 X commute (but not and d 3 X !) we obtain: 



dtVpjt) _ 
Vv(t) V v (t) 



(12d) 

Introducing the scaled (t-)expansion 9 := NO we define an effective (t-)Hubble function 
in the hypersurfaces by 

(«)„ = = S*2H =: Mb . (12e) 

Vt> {t) Civ 

(Notice that we reserve the overdot for the covariant derivative.) 
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It is now straightforward to prove the following Lemma for an arbitrary scalar field T(X l , t): 
Lemma (Commutation rule) 



d t (r)v - (d t r) v = (re)v - (y)v(e) v , (12/) 

or, alternatively, 

d t (r) v + 3H v (y) v = (d t r + re) v . (i2 9 ) 

A simple application of this Lemma is the proof that the total restmass in a domain 
is conserved: let T = g, then dt{g)v + 3Ht>(q)t> = (dtQ + q9)t> = according to the local 
conservation law Eq. (9e). q.e.d. 

3.2. Averaged Equations for Irrotational Perfect Fluids 

Averaging Raychaudhuri's equation (10b) and the Hamiltonian constraint (7a) with the 
help of the prescribed procedure, we end up with the following two equations for the scale 
factor of the averaging domain which we may formulate in the form of a theorem: 

Theorem — Part I (Equations for the effective scale factor) 

The spatially averaged equations for the scale factor ax>, respecting restmass conservation, 
read: 

averaged Raychaudhuri equation for the scaled (t-) densities e := N 2 e and p := N 2 p: 

3^E + 4 nG (e + 3p) v = Q v +Vv ; (13a) 
ax> 

averaged Hamiltonian constraint: 

ml - 167rG(£)x, = -(Qv + (K)t>) , (136) 

where we have introduced the scaled spatial (t-)Ricci scalar 1Z := N 2 1Z, and we have 
separated off the domain dependent 'backreaction' terms: the kinematical backreaction , 

Qv:=[2(N 2 ll) v - 2 -(Nl)l\ = |((gf - (0) v )\ - 2(a 2 ) v , (13c) 

with the scaled (t-) shear scalar a := Ncr, and the dynamical backreaction , 

Vv := (A)v + {N9)v , (13d) 
with the scaled (t-) acceleration divergence A := N 2 A. 
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Note: Eq. (13a) can also be obtained by an argument given by Yodzis (1974), which is 
summarized in Appendix C of (Paper I). 

The source terms on the r.-h.-s. of Eq. (13a) include the 'kinematical backreaction' (13c) 
that describes the impact of inhomogeneities on the scale factor due to averaged shear 
and expansion fluctuations. It vanishes for the standard FLRW cosmologies. Additionally, 
Eq. (13a) features another 'dynamical backreaction' term (A)t> together with a technical 
term due to the change of the lapse function. The former term also vanishes for stan- 
dard FLRW cosmologies; both terms vanish for zero pressure. Note that the averaged 
Hamiltonian constraint does not involve pressure terms as expected. 

These equations show that the averaged shear fluctuations tend to increase the expansion 
rate similar to the effect of the averaged energy source terms (provided the energy condition 
(e + 3p)x> > holds), while the averaged expansion fluctuations work in the direction of 
stabilizing structures. Pressure forces can do both; the sign of the averaged divergence of 
the 4-acceleration can be positive or negative. In the Newtonian framework one can show 
that, to a first-order approximation, the combined effect of gravity and pressure leads 
to stabilization of structures (Buchert & Dommguez 1998, Buchert et al. 1999, Adler & 
Buchert 1999). Note, however, that since pressure is a source of the gravitational field 
energy too, it is harder to oppose the gravitational collapse than in the corresponding 
Newtonian treatment (compare the terms which add positive contributions in Eq. (9d) 
with their Newtonian analogues). 

We proceed by calculating the integrability condition for the system of equations (13a,b), 
i.e., we shall answer the question which equation has to hold in order that (13b) be the 
integral of (13a). For this end we take the partial time-derivative of (13b) and insert into 
the result again our starting set of equations (13a) and (13b). We get: 

Theorem — Part II (Integrability and energy balance conditions) 
Eq. (13b) is an integral of Eq. (13a), iff 

d t Qv + QH V Q V + dt(TZ)v + 2H V (K) V + AH V V V 

-16TrG[d t (e)-D + 3H-d(s + P)-d] = . (14a) 

The expression involving the energy density and the pressure does not vanish in general. 
To see this we average the local energy conservation law (9a) . We obtain: 

d t (e)< D + 3H< D (e + p)< D = (d tP )< D -d t (p)v • (146) 
For the scaled (t-)variables we accordingly have for the local law: 

d t e + e(e + p) = 2Ne , (14c) 

and for the average: 

d t (e} v + 3H v (e + p} v = (d t p} v -d t (p) v + (2Ne} v . (Ud) 
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This shows that the pressure term introduces a possibly interesting effect. In the 'dust' 
case (Paper I: Corollary 1) the averaged fields obey the same equations as the local fields 
provided we use their representation in terms of invariants of the second fundamental 
form. Here, this is no longer true. In particular, Eq. (14b) shows that the averaged energy 
conservation law invokes non-commuting terms that are nonzero for inhomogeneous fluids. 
Thus, even if both the 'kinematical' and 'dynamical backreaction' terms are assumed to be 
negligible or cancel for "some" reason, the averaged model is different from the standard 
homogeneous-isotropic models. This fact will be manifest in the following more compact 
alternative representations of the averaged equations. 



Corollary 1 (Averaged equations: first effective form) 
Let us define effective densities as follows: 

Then, the averaged equations can be cast into a form similar to the standard Friedmann 
equations: 

3^ + 4^(^ + 3^) = ; (15c) 

ml + - 16ttGs^ = , (15d) 

and the integrability condition of (15c) to yield (15d) has the form of a balance equation 
between the effective sources and the averaged spatial (t-)Ricci scalar: 

d t e^ + 3H v {e^+p^) = ^(^(7^ + 2^(7^) . (15e) 

The effective densities obey a conservation law, if the domains' curvature evolves like in 
a "small" FLRW cosmology, (JZ)t> = 0, or (JZ)v oc a^, 2 , respectively. In particular, Eq. 
(15e) shows that in general the averaged densities are directly coupled to the evolution of 
the averaged spatial curvature. 

Considering the averaged spatial t-Ricci scalar as an effective source as well, one may 
cast the equations into an even more elegant form. 
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Corollary 2 (Averaged equations: second effective form) 



Defining 

(2) ,~ x Q.V {K)t> 

£ °* :={£) » - - IfoFG ' ^ (16a) 

(2) Qd (7^)x» Pr> , 1ft ,x 

^^==^"16^ + 48^-12^ ' (166) 

we obtain equations that assume the form of spatially 3-Ricci flat Friedmann cosmologies: 



ax> 



+ AttG (4g + 3p£>) = ; (16c) 



6#|, - lGirGe^ = , (16d) 

and the integrability condition of (16c) to yield (16d) has exactly the form of a conservation 
law: 

cVg+Stfx^+rf?) = • (16c) 

Remarks: 



These alternative representations reduce the solution of the averaging problem for scalars, 
at least formally, to the problem of finding an 'effective equation of state' that relates the 
effective densities. Relativistic Lagrangian perturbation schemes will be useful to establish 
such relations. Looking at Eqs. (15) it is interesting to note that the 'kinematical backre- 
action' term itself effectively performs like a free scalar field source, or like 'stiff matter' 
in the case (TZ)t> oc a^, 2 . However, care must be taken with such statements, since Qx> is 
related to (TZ)x>, and there is no a priori reason why {TZ)x> oc a^ 2 if backreaction is present. 
It nevertheless suggests to separate off the 'stiff component' from the 'effective equation 
of state': we already noted that — 16 * G Qt> in Corollary 1 forms a 'stiff' part; deviations 
from 'stiffness' are, apart from those due to the matter sources, due to the_ 'dynamical 
backreaction' — ^q Vt>- In the form of Corollary 2 the 'stiff' part is — 16 ^ g (2t> + (fi)v)-> 
and the deviations from 'stiffness' are due to the term j^q((T^)v — *Pt>)- If 'dynamical 
backreaction' compensates the averaged scalar curvature, then the whole 'backreaction' 
forms a 'stiff component'. The condition for this can be inferred from Eqs. (13a,b) yielding 
the general relation: 

(n}v-V v = -3(^ + 2Hl) + 12nG(e-p} v . (17) 



In general, the system of equations (13,14) is not a closed system, which can be most 
easily seen in the form of Corollary 2: we have three equations (16c, d,e) for the three 

(2) (2) 

variables a-p, e e ^ and p eii ', but only two of them are independent. We need an effective 
equation of state to close the system. 
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4. Discussion of Subcases 



Since this paper is meant to provide the basic architecture for applications, let us note the 
following useful formulas. 

Firstly, the equations of Sect. 3 simplify by using the following reparametrization of 
time: the line element is invariant under the change of the time coordinate t h- > S(t), so 
there is still some gauge freedom. Using the 'phase fronts' as the new time coordinate we 
define a new lapse function N by 

Ndt=:NdS , i.e. , N = Nd t S . (18a) 

In particular, the total (Lagrangian) derivative becomes 

where the latter equality follows from Eqs. (lc) and (4b). 

Notice that with this new choice of time coordinate all equations of Sect. 3 remain 
form invariant, if iV is replaced by N, and partial time-derivatives are replaced by partial 
derivatives with respect to S. The latter will be abbreviated by a prime in what follows. 
All fields are functions of the independent variables (X 1 , S) now. 

Secondly, for the evaluation of the terms appearing in the averaged equations it is 
helpful to note the following simplifications. We shall give expressions for a barotropic 
fluid and, especially, for the simple class of matter models p = je, which is relevant for 
many applications. 

For barotropic fluids, 

e + 9(e + a(e)) = , (18c) 

we can integrate the energy conservation law along trajectories of fluid elements using 
J = 9J (Eq. 12c) and g = —9g (Eq. 9e) to find the entropy density, 

/de 
-— , (18c?) 
e + a[s) 

and, upon performing the integral, the energy density. In particular, for a' = 7 = const. 
we obtain: 

s(e)cce^ ; e oc . (18e) 

Hence, with g oc J -1 , and using Eq. (llg), 

e oc g 1+1 ; h oc g 1 ; s oc g . (18/) 

Furthermore, the scaled (t-)variables are now normalized with respect to the magni- 
tude h or its square, respectively. E.g., we have 

N=l ; ~e = N6= 6 - ; e=^ ; etc. . (18a) 
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For barotropic fluids with a' = 7 = const, the normalization function can be written in 
powers of the restmass density 3 . Also for this case the expression involving the covariant 
derivative of the lapse function in Eq. (13d) is simply proportional to the t— expansion rate 
(upon using the integral for h Eq. (18f)): 

~ N' h' o' 

n = -L = ~ = - 1 ^- = 1 e . (i8/i) 

N h q 

The term involving the change of the lapse function in Eq. (14d) is simply time-dependent 
in two cases: it can be written as follows for the homogeneous case (Eq. (18i); Subsect. 4.1), 
and for the 'stiff fluid' representing a free minimally coupled scalar field (Eq. (18j); Sub- 
sect. 4.4): 

(2Ne} v = GH<J^- , (18i) 

(29e)v = 3Hv ■ (18j) 

It is interesting that the latter (fully inhomogeneous) case implies simplifications such that, 
e.g., Eqs. (14c, d) reduce to identities (since in this case e = p = |, compare Subsect. 4.3). 

For the 'dynamical backreaction' we have for p = 75: 

Vv = (A}v + l(e 2 }v , (18k) 

and, noting that 6 = — and using Eqs. (10a) and (18f) with h = C1Q 1 , the 'dynamical 
backreaction' term can be entirely written in terms of the restmass density: 

The same is true for the averaged spatial t— Ricci curvature using Eqs. (10c) and (18f) 
with e = C 2 Q 1+1 : 

11 

(K) v = 12wG^(1- 1 ){q 1 -^) v -q(^J +V v , (18m) 

where from Eq. (lie) C\ = (1 + 7)6*2, and C 2 is determined by initial conditions. This 
makes the problem accessible to relativistic Lagrangian perturbation models, since the 
restmass density can be integrated exactly along the flow lines, g oc J -1 , and J can be 
computed from the basic dynamical variable in Lagrangian perturbation theory (see, e.g., 
Kasai 1995, Takada & Futamase 1999 for relativistic 'dust' models, Adler & Buchert 1999 
for pressure-supported fluids in Newtonian theory). 



3 The constants C\ and C 2 appearing in the following equations are in general X 1 - 
dependent; also equations of state, if they arise as integrals of dynamical equations, involve 
X*-dependent functions of integration such as 7. One may use the freedom to relabel the 
fluid elements such that the constants or some product of them are ^'-independent. We 
here assume that all the constants are equal for each fluid element for notational ease; it 
is straightforward to write down the more general expressions, if needed. 
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4.1. Homogeneous— isotropic Cosmologies 



The requirement of homogeneity and isotropy reduces Eqs. (13a,b) to the familiar Fried- 
mann equations. We may simply put the lapse function (without loss of generality) equal 
to 1 in Eqs. (13a,b) and find Q.v = Vt> = 0, and: 

3- +4tvG(e h + 3p H ) = ; (19a) 
a 

6H 2 + TZ H -16nGe H = ; H := - . (196) 

a 

Already before averaging eu and pu are functions of time only, and the scale factor assumes 
its global standard value a-p = a. The domain dependence has disappeared. Note also that 
the integrability condition (14a) reduces to the equation for the 3-Ricci curvature scalar 
of the spatial hypersurfaces, 

7?o 

n H + 2Hn H = o n H = -^- , (i9c) 

a z 

and the averaged conservation law (14b) coincides with the local one, 

e H + 3H(e H +p H ) = . (19d) 

For a given relation between pu and en the system of equations (19) is closed. 

Alternatively, we may look at the homogeneous-isotropic models within the present 
framework in terms of a time-dependent lapse function. For the time variable S and 
requiring h = hn(S), e = en(S), p = pn(S), 1Z = TZh(S), the domain dependence 
disappears, and we have N = hjj 1 , 9 = 0H(S)hJ I 2 1 9h = 3H. Setting a = we also have 

Qd = 0, but Vd = —j^Oh ^ 0. The system of averaged equations Eqs. (13) together 
with Eqs. (14) reduces to the following set: 

3— + 4ttG^ t (s h + 3ph) = -3H^ l , (20a) 
a h H hu 



H / H a H n H 

g H + 3Hq h = , (20d) 



which, together with hn = £H ^ H an d an equation of state pn = ct(en), are four equations 
for the four unknown functions a(S), Qh(S), Eh(S), and TZh(S), but only three equations 
are independent. We have to use the fact that, with pu = a(en), £h can be expressed in 
terms of qh, which closes the system. 
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This set of equations looks odd with respect to the hypersurfaces S = const. However, 
upon reinvoking the covariant time derivative, e.g., d = /i#a', H = ^ /i^ 1 , hu in the first 
three equations disappears, and we recover the familiar form of Eqs. (19), illustrating the 
covariance of the averaged equations. 

4.2. Inhomogeneous 'Dust' Cosmologies 

Putting in Eqs. (13) all pressure terms to zero and noticing that e reduces to the restmass 
density we have h = 1, the lapse function N = 1, and the covariant time-derivative 
= ds- Hence, we directly recover the form of the averaged equations of Paper I for 
cosmologies with a 'dust' matter content: 

3— +4irG(g} v = Q v ; (21a) 

Clj) 

6Hl + (JI) v -16itG(q)v = -Qv ; H := — , (216) 

(It* 

with the integrability condition Eq. (14a) being 

Qb + QH v Q v + {TZ)h + 2H v {TZ) v = 16tvG((q)v + 3Ht>(q)v) ; (21c) 

the balance equation (14d) reduces to the continuity equation for the averaged restmass 
density: 

(Q)h + ^H v ( e ) v = . (21d) 

Notice that only in the case of an inhomogeneous 'dust' model we can put the lapse 
function N or N itself equal to 1 without loss of generality; the averaged equations are 
already covariant in the comoving and synchronous gauges. It is to be emphasized that 
the 'dust matter model' cannot generically be foliated into hypersurfaces S = const, with 
an inhomogeneous lapse, for Eqs. (9) necessarily imply a constant lapse function for the 
geodesic condition of vanishing acceleration that itself is implied by vanishing pressure. 

Eqs. (21) form a set of four equations for the four unknown functions ap, (q)t>, 
(1Z)t>, and Qx>, but only three equations are independent. As discussed in Paper I, this 
system cannot be closed on the level of ordinary differential equations unless additional 
(e.g. topological) constraints are imposed. Of course, this remark also applies to the more 
general matter models. 

From Eqs. (21) it is obvious that the requirement Qx> = is necessary and sufficient 
in order that ax>(t) obeys the equations of standard FLRW cosmologies. 



16 



4.3. Radiation— dominated Inhomogeneous Cosmologies 



Let us consider a situation in which radiation is directly coupled to the matter fluid (ac- 
cording to the conjecture of a local thermodynamic equilibrium state of radiation and 
matter) , then we may describe the radiation cosmos as a single component perfect fluid 
with radiation pressure r p 1 and radiation energy density £ 7 obeying e 7 = 3p 7 (Ellis 1971, 
see, however, Ehlers 1971). We infer already from the averaged conservation law (14b) 
that the time evolution of a radiation-dominated inhomogeneous universe is different from 
that expected from the corresponding homogeneous-isotropic world model: 

d t (ej)v +AH v {e 1 ) v = (dtP 7 )x> - &(p 7 )x> ■ 

The term on the r.-h.-s. of this equation, which vanishes in the standard model, may be 
interpreted as an accumulated effect from inhomogeneities in the radiation field yielding 
deviations from a global 'equilibrium equation of state'. It should be stressed that these 
deviations also occur in the case where the 'backreaction' terms Qt> and Vd are both 
found or assumed to be negligible. Therefore, radiation-dominated fluids deserve further 
detailed study. 



4.4. Inhomogeneous Scalar Field Cosmologies 

Following Bruni et al. (1992) we may describe the dynamics of a scalar field 0, mini- 
mally coupled to gravity, in terms of the natural slicing of spacetime into a foliation of 
4> = const, hypersurfaces. Einstein's equations for a scalar field source are (under condi- 
tions stated below) equivalent to the phenomenological 3 + 1-description of an evolving 
pressure-supported perfect fluid with energy-momentum tensor and corresponding fluid 
4-velocity (normal to the hypersurfaces of constant 0) (Taub 1973, Madsen 1988): 

T /t = £<t>UnU v + p^h^ ; ulJl = ■ ( 22a ) 

The magnitude ip normalizes the momentum density vector <9 M so that u^u^ = —1, 

< V = y/- da <f> d <*<f> = u^d^ =: $ , (226) 

where the overdot stands for the material (or Lagrangian) derivative operator as before. 
From Eq. (22a) we conclude that Einstein's equations feature the perfect fluid energy- 
momentum tensor (e.g., Madsen 1988): 

T U = ( £ <f> + Pd>)u^u u + p^g^ u = d^(f)d v (f) - Qd a 0<9 a + V eS ((j))^ , (23a) 



with 



e^ = ^ 2 + V eS (<P) , (236) 



17 



V4> = 2^ ~ VeS ^ 



(23c) 



The fluid analogy is valid, if the fluid 4-velocity is time-like and, hence, the hypersurfaces 
4> = const, space-like. For this to be true the 4-gradient of the scalar field has to be 
time-like, 

-V> 2 = d a (f)d a (t) < . (23d) 
This is a sufficient requirement for having the energy condition 

+ p^ = ^ 2 > , (23e) 

which follows from Eqs. (23b, c). This condition still allows for powerlaw inflation; ex- 
ponential inflation is excluded and has to be studied as a separate case. This case can 
be studied within the fluid analogy, if we model the constant effective potential with a 
cosmological constant. (The basic equations have to be used including the cosmological 
constant - see Appendix). 

We finally note that, only for constant effective potential (V^ ff = 0), the restmass 
conservation law of a perfect fluid corresponds to the Klein-Gordon equation 

V, + 0V> + Kff(0) = O , V> = • (23/) 

The 'equation of state' of a scalar field is, in general, not barotropic (see: Bruni et al. 
1992). However, for interesting cases it is barotropic and can be represented in terms of an 
'equation of state' p^ = a ( f > (e ( f > ): 1) = —1 for a stationary state (vacuum ground state), 
and 2) = +1 for the free state (corresponding to a 'stiff fluid'). In general, if there 
exists an equation of state, it will have the form p^ = a^e^, s^) with the entropy density 
Sff,. However, an evolving scalar field will in general yield a dependence of p^ on the other 
dynamical variables and 0. The function a<j, is determined by the dynamics and it may 
or may not be a priori written, e.g., as a function of the density and the entropy density. 

Since the minimally coupled free scalar field (dilaton) is singled out in the present inves- 
tigation as the only inhomogeneous case in which the averaged equations attain a simpler 
form, it is certainly worth studying this case in more detail. This is the headline of a 
forthcoming work (Buchert & Veneziano 2001). 
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Appendix: 

Basic Equations in the ADM Formalism 

Let be the future directed unit normal to the hypersurface E. The projector into E, 
= g^ u + n^n v , (=>• h^n^ = , h} i v h v 1 = h^), induces in E the 3-metric 

hij := g^h^j . (Ala) 

Let us write 

n M = iV(-l, 0,0,0) , ^ = ^(1,-AH , (Alb) 

with the lapse function N and the shift vector N\ Note that N and N l can be determined 
by the choice of coordinates. 

From n M = g^ v n v we find #00 = -(N 2 - iVjAP); goi = Nf, gij = h t j and, setting x° = t, 
the line element becomes: 

ds 2 = -(N 2 -N l N l )dt 2 + 2N l dtdx l +g lj dx i dx j = -N 2 dt 2 +g lJ (dx l + N l dt)(dx J dt) . 

(Ale) 

Introducing the extrinsic curvature tensor on E by 

Kij := -n il . v h ii i h v j = -n t .j , (Aid) 
we obtain the ADM equations (Arnowitt et al. 1962, York 1979): 

Energy (Hamiltonian) constraint: 

71 - K l 3 K\ + K 2 = 16nGe + 2A , e := T^n u ; (A2a) 

Momentum constraints: 

K* mi -K\\ 3 =MGJ 3 , Ji:=-T^h\ ; (A2b) 

Evolution equation for the first fundamental form: 

jjdtgij = -2Kij + l(JV iNj . + N jni ) ; (A2c) 
Evolution equation for the second fundamental form: 

IfttfV = TV, + KK*, - S^A - ±N\% + 1 (k\N% - K k jN \ k + N k K^ k ) 

-torG&j + U l j(e - S k k )) , (A2d) 
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where S ZJ := T^h^h^ . 

For the trace parts of (A2c) and (A2d) we have: 

±d t g = 2g{-K+±N\ k ) , g := detfoy) ; (A2e) 
±d t K = n + K 2 -4nG(3e-S k k )-3A-±N\\ k llk + ±N k K llk . (A2f) 

For our purpose of averaging we have used equations that correspond to the coordinate 
choice of vanishing shift vector. Thus, all inhomogeneities of the fluid were put into the 
3-metric and the lapse function. 

Assuming the tensor T^ v has the form = eu^u v + ph^ v and putting the shift vector 
N l = and also A = 0, we obtain the equations of the main text by defining the lapse 
function in such a way that a« = = -j^j- = — -jp. 

Notice that with this choice the unit normal coincides with the 4-velocity and, especially, 
the momentum flux density in E vanishes. The total time-derivative operator of a tensor 
field T along integral curves of the unit normal, jjrJF := n v d v T = u v d v T becomes 

since n v T\\ v = 0. Note that, although the definition of proper time is r := / Ndt, the line 
element cannot be written in the form of the comoving gauge by measuring "time" through 
proper time dr = Ndt, since dr is not an exact form in the case of an inhomogeneous lapse 
function. The exterior derivative of the proper time will involve a non-vanishing shift 
vector according to the space-dependence of the lapse function. Therefore, a foliation into 
hypersurfaces r = const, with simultaneously requiring u a = —d a T is not possible. 

For vanishing shift vector the line element reads: 

ds 2 = -N 2 dt 2 + g ij dX i dX j . (A3b) 
The lapse function itself may be written explicitly: from h = S we have: 

N = \d t S . (ASc) 

Note that, if we assume that h > 0, implying the energy condition e + p > 0, the proper 
time advances only in periods when the derivative dfS > which makes a difference if we 
consider fluids that mimick a scalar field source. 

The coordinates in E are written in capital letters now, because for vanishing shift vec- 
tor they correspond to Lagrangian coordinates as in classical fluid mechanics. In these 
coordinates £ = ±f t = jfd t \ x *- 
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